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The ab initio no-core shell model (NCSM) is extended to include a realistic three-body interac- 
tion in calculations for p-shell nuclei. The NCSM formalism is reviewed and new features needed 
in calculations with three-body forces are discussed in detail. We present results of first applica- 
tions to 6 ' 7 Li, 6 He, 7 ' 8 ' 10 Be, 10 > n ' 12 B, 12 N and 10 .ii,i2,i3 C uging the Argonne V8' nucleon-nucleon 
(NN) potential and the Tucson-Melbourne TM'(99) three-nucleon interaction (TNI). In addition 
to increasing the total binding energy, we observe a substantial sensitivity in the low-lying spectra 
to the presence of the realistic three-body force and an overall improvement in the level-ordering 
and level-spacing in comparison to experiment. The greatest sensitivity occurs for states where the 
spin-orbit interaction strength is known to play a role. In particular, with the TNI we obtain the 
correct ground-state spin for 10 > n B and 12 N, 12 B, contrary to calculations with NN potentials only. 

PACS numbers: 21.60.Cs, 21.45.+V, 21.30.-x, 21.30.Fe 

I. INTRODUCTION 

The ab initio no-core shell model (NCSM) Q is a method to solve the nuclear structure problem for light nuclei 
considered as systems of A nucleons interacting by realistic inter-nucleon forces. The calculations are performed using 
a large but finite harmonic-oscillator (HO) basis. Due to the basis truncation, it is necessary to derive an effective 
interaction from the underlying inter-nucleon interaction that is appropriate for the basis employed. The effective 
interaction contains, in general, up to ^4-body components even if the underlying interaction had, e.g. only two-body 
terms. In practice, the effective interaction is derived in a sub-cluster approximation retaining just two- or three-body 
terms. A crucial feature of the method is its convergence to exact solution with increasing basis size and/or an increase 
in the effective interaction clustering. 

In the past, applications were limited to only realistic two-nucleon interactions. However, the recent introduction of 
the capability to derive a three-body effective interaction Q and apply it in either relative-coordinate Q or Cartesian- 
coordinate 4] formalism together with the ability to solve a three-nucleon system with a genuine three-nucleon force 
in the NCSM approach 5] opens the possibility to include a realistic three-nucleon interaction (TNI) in the NCSM 
Hamiltonian and perform calculations for the p-shell nuclei. 

We note that there are several methods that can be used to solve the A = 3, 4 systems with realistic Hamiltonians 
that include a realistic three-body interaction. However, until now, only the Green's function Monte Carlo (GFMC) 
method [a, [3, la, |9j is capable to obtain solution for light p-shell nuclei with a Hamiltonian that includes both realistic 
two- and three-nucleon force. 

In this paper, we introduce an extension of the NCSM formalism to accommodate the TNI and present first 
applications for several p-shell nuclei. The main purpose of this paper is to introduce the formalism needed to 
perform ab initio shell-model calculations for the p-shell nuclei with Hamiltonians that include a realistic three- 
nucleon interaction. At the same time, we present a snapshot of first applications for several p-shell nuclei of different 
masses with the primary goal of assessing a general impact of a realistic TNI on the structure of different p-shell 
nuclei. In this study, we limit ourselves to the use of a single TNI, the chiral-symmetry based Tucson-Melbourne 
TM'(99) [13, combined with the Argonne V8' NN potential @. More detailed studies using a broader variety of 
realistic three-nucleon interactions will follow in the future. 

In Sectional the extension of the NCSM formalism to include realistic three-nucleon forces is discussed. In Section 
IIIII we present first results obtained with Hamiltonians that include the Tucson-Melbourne TM'(99) TNI for 6:7 Li, 6 Hc, 
7 - 8 ' 10 Be, io.h.12b ) i2 N an( j 10,11,12, 13q g ome overa n observations of the effect of the TNI are gathered in Section llVl 
Finally, we summarize our conclusions in Section Ivl 

II. AB INITIO NO-CORE SHELL MODEL WITH A THREE-NUCLEON FORCE 

A detailed description of the NCSM approach was presented, e.g. in Refs. Here, we emphasize extensions 

and modifications needed when a genuine TNI is included. In the case when the TNI is considered, the starting 
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where m is the nucleon mass, T^JN.y , the NN interaction with both strong and electromagnetic components, VNNN,ijfc 
the three-nucleon interaction. In the NCSM, we employ a large but finite harmonic-oscillator (HO) basis. Due to 
properties of the realistic nuclear interaction in Eq. , we must derive an effective interaction appropriate for the 
basis truncation. To facilitate the derivation of the effective interaction, we modify the Hamiltonian JQ) by adding to 

it the center-of-mass (CM) HO Hamiltonian H CM = T C m + U C m, where U C m = \AmVL 2 R 2 , R = \ Y,i=\ The effect 
of the HO CM Hamiltonian will later be subtracted out in the final many-body calculation. Due to the translational 
invariance of the Hamiltonian Q the HO CM Hamiltonian has in fact no effect on the intrinsic properties of the 
system in the infinite basis space. The modified Hamiltonian can be cast into the form 
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Next we divide the ^4-nucleon infinite HO basis space into the finite active space (P) comprising of all states of up to 
iV max HO excitations above the unperturbed ground state and the excluded space (Q = 1 — P). The basic idea of the 
NCSM approach is to apply a unitary transformation on the Hamiltonian J2J), e^ s H^e s such that Qe~ s H^e s P = 0. 
If such a transformation is found, the effective Hamiltonian that exactly reproduces a subset of eigenstates of the 
full space Hamiltonian is given by H c g — Pe~ s H^e s P. This effective Hamiltonian contains up to A-body terms 
and to construct it is essentially as difficult as to solve the full problem. Therefore, we apply this basic idea on a 
sub-cluster level. When a genuine TNI is considered, the simplest approximation is a three-body effective interaction 
approximation. The NCSM calculation is then performed with the following four steps: 

(i) We solve a three-nucleon system for all possible three-nucleon channels characterized by a fixed angular momen- 
tum, isospin and parity with the Hamiltonian J5Jl, i.e., using 



Ti.3 = hi+h 2 + h 3 
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Note that in Eq. ©, A > 3 corresponding to the investigated nucleus is kept in V^' A ,i,j = 123 as defined m 
© ■ Consequently, the three nucleons feel a pseudo-mean field of the spectator nuclcons generated by the HO CM 
potential. It is necessary to separate the three-body effective interaction contributions from the TNI and from the 
two-nucleon interaction. Therefore, we need to find three-nucleon solutions for the Hamiltonian with and without the 
Vnnn.123 TNI term. The three-nucleon solutions are obtained by procedures described in Refs. 0,0 (without TNI) 
and ja] (with TNI). In particular, we note that we use two-body effective interaction corresponding to a very large 

(3) 

three-nucleon basis space, i.e., up to A^ax = 40, to obtain these solutions. The number of three-nucleon channels to 
include in the calculations is restricted by the maximal active space N max to be used in the final A-nucleon calculation. 

(ii) We construct the unitary transformation corresponding to the choice of the active basis space P3 from the three- 
nucleon solutions using the Lee-Suzuki procedure ^30, El EH- I n particular, we demand that Q3e~ s Hze s P3 = 0. 
Let us remark that the definition of the three-nucleon active space P3 is determined by the definition of the A-nucleon 
active space P. The P3 states contains all three-nucleon states up to the highest possible three-nucleon excitation, 
which can be found in the P space of the A-nucleon system. For example, for A = 6 and iV max = 6 (6h~il) space 
we have P3 defined by ./V3 max = 8. Similarly, for the p-shell nuclei with A > 7 and -/V max; = 6 (6hQ) space we have 
^3max = 9. We note that the anti-hermitian transformation operator S can be expressed as | l 1| 



S = arctanh(o; — ur ) 
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with the operator u> satisfying u> — Q3OJP3. If the eigensystem of the Hamiltonian H3 © is given by 7Y3 1 fc) = Ek\k), 
then the operator uj can be determined from 



(a Q \uj\ap) = ^2(a Q \k)(k\a P ) 
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where the tilde denotes the inverse of the matrix defined by matrix elements (ap\k), i.e., ^2 ap (k\ap }(ap\k / ) = 5k,k' 
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for k, k' E JC. In Eq. JSJ, \ap) and \o.q) are the active-space and the Q-space basis 



3 



states, respectively, and K, denotes a set of dp eigenstates, whose properties are reproduced in the P-space, with dp 
equal to the dimension of the P-space. The set K. is typically chosen by selecting the lowest three-nucleon eigenstates 
for the given three-nucleon channel. In practice, to calculate the P-space matrix elements of the effective Hamiltonian, 
we apply the Eqs. (15) and (16) of Ref. Q. The three-body effective interaction is then obtained as 
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Note that we distinguish the transformation Snn+nnn corresponding to the Hamiltonian TL 3 IjSJ and the transforma- 
tion Snn corresponding to the Hamiltonian without the TNI term, i.e., hi + h% + h 3 + V 12 ' + V 13 + V 23 . The 
three-body effective interaction contribution from the TNI we then define as 

t/NNN _ T/ NN+NNN _ x/NN / s \ 
|/ 3eff,123 — V 3eff,123 *3eff,123 • \° I 

(iii) As the three-body effective interactions are derived in the Jacobi-coordinatc HO basis but the p-shell nuclei 
calculations will be performed more efficiently in a Cartesian-coordinate single-particle Slater-determinant M-scheme 
basis, we need to perform a suitable transformation of the interactions. This transformation is a generalization of the 
well-known transformation on the two-body level that depend on HO Brody-Moshinsky brackets. The antisymmetrized 
Jacobi-coordinate HO basis in which the three-body effective interaction is obtained can be represented in the form 
\NiJMTMx) where N is the total number of HO excitation, i is an additional quantum number enumerating the 
antisymmetrized states, J and T are the total three-nucleon angular momentum and isospin, respectively, and M, 
Mt their third components. This basis can be expanded in a 2 + 1 nuclcon cluster basis using expansion coefficients, 
i.e., 

\NiJMTM T ) = ^((nZs^,^£^ l /i)||iViJT)|(ri/sit,7V'£i l 7i)JMTM T ) , (9) 

where the \nl) HO state describes relative motion of nucleon 1 and 2, \AfC) describes the relative motion of the third 
nucleon with respect to the CM of nucleons 1 and 2, and ((nlsjt.N C^J^)\\NiJT) is the coefficient of fractional 
parentage. It holds that N = 2n + I + 27V + C. For further details and the procedure how to calculate the expansion 
coefficients sec Ref. 3J. The basis for the three-body effective interaction suitable for the shell- model code input is 
a three-nucleon Slater-determinant HO basis characterized by M 3 = mj 1 + nij 2 + m 33 , Mt = mt 1 + rn t2 + m t3 and 
parity, formed by the single particle states \nl^jmj^mt) that depend on the single-particle coordinates. We introduce 
a short-hand notation for these Slater-determinant states, i.e., K^-gi^- Am t )( a (, c )), where a, b, c labels the occupied 
states. In order to transform the three-body effective interaction to the new basis we need first to couple the relative 
coordinate basis JjJJ with the three-nucleon CM HO states |7VcmPcm7Wcm) to form a complete basis. The three-body 
effective interaction is independent on the A^cm-^cm-^cm quantum numbers. The overlap of the two respective states 
can then be expressed in the form 

{(nl-jmj-m t )( abc )\NiJMTMT; N C mL C mM C m) = fa 
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where we explicitly show the total HO quantum number conservation that also implies parity conservation, the angular 
momentum and isospin third component conservation. The sum goes over the expansion A, A, A12, L12 and the m 
quantum numbers in the Clebsch-Gordan coefficients not appearing on the left-hand side. The (nilin2l2L\n3h n 4hL)d 
are the generalized Brody-Moshinsky brackets for two-particles of mass ratio d as defined e.g., in Ref. 15]. Sums of 
the Clebsch-Gordan coefficients in Eq. 1)11) [) can be re-expressed as sums of 6-j and 9-j coefficients. However, some 
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Clebsch-Gordan coefficients will remain as we have a magnetic quantum number dependence on the left-hand side. It 
should be pointed out that the overlap (|10|l is independent of the HO frequency f2, and for a given N maX) it needs to 
be calculated just once. The three-body effective interaction M-scheme matrix elements are then obtained by 

((71/ ^cs^^nl-jriij-m^tdef)) = ^2 (( nl 2~j m j2 m t)(abc)\NiJMTMT]NcMLcMMcM} 

(NiJT\V 3cSA23 \N'i'JT) {N'i'JMTM T -N CM L CM M C M\{nl^m^rn t ) {de}) ) . (11) 
(iv) We solve the Schrodinger equation for the A-nucleon system using the Hamiltonian 

A ^ A A 

H A,cS = E ^ + A _ 2 E V teG,ijk + E V 3cS%k > ( 12 ) 
i—1 i<j<k i<j<k 

where the -^z^ factor takes care of over-counting the contribution from the two-nucleon interaction. At this point 

we also subtract the Hqm and add the Lawson projection term P(Hqm — ffifi) to shift the spurious CM excitations. 
Just as the bare interaction is translationally invariant, so is the effective interaction, and the energies of the physical 
eigenstates corresponding to the Ofi.fi excitation of the CM are independent on the choice of (3. The A-nucleon 
calculation is then performed using either the Many-Fermion Dynamics shell model code |l(l| generalized to handle 
three-body interactions or using the newly developed code REDSTICK [Tj| , using a different algorithm for evaluating 
the Hamiltonian matrix elements, with the Hamiltonian in the form: 

A A 

+ E v ™«" + E v w + - \m . (is) 

i<j<k i<j<k 

We summarize this section by repeating that the effective interaction depends on the nucleon number A, the 
HO frequency f2, and the P-space basis size defined by A^ max . It is translationally invariant and with A max — > oo 
the NCSM effective Hamiltonian approaches the starting bare Hamiltonian (|TJ. Consequently, with the basis size 
increase the NCSM results become less and less dependent on the HO frequency and converge to the exact solution. 
Alternatively, by increasing the clustering of the effective interaction for a fixed P-space size the NCSM effective 
Hamiltonian approaches the exact ^4-nucleon effective Hamiltonian that reproduces exactly a subset of eigenstates of 
the starting Hamiltonian JIJ. 



H 



o,sm 

A,eS 



E 



(Pi ~Pj) mQ 



2Am 



III. APPLICATION TO SELECTED p-SHELL NUCLEI 



The calculations with the three-body interactions are technically and computationally intensive. When a genuine 
TNI is employed the step (i) from the previous section is computationally demanding. In addition, when a three-body 
effective interaction is used in general, the steps (iii) and (iv) are complex as well. In these first NCSM calculations 
with the TNI in the p-shell, we limit ourselves to the use of a single realistic interaction. In particular, we use the 
isospin invariant AV8' NN potential |(| that is a simplified version of the AV18 NN potential. For the TNI, we employ 
a version of the original Tucson-Melbourne force 18] corrected to drop a term later found to be inconsistent with 
the chiral-symmetry-based effective field theory TNI of two-pion range [19j. We note that there have been different 
strength constants (taken from pion- nucleon scattering data) suggested for this force 0, 0, 0] . We use only one 
particular set, i.e. the TM'(99) updated strength constants of the Tucson-Melbourne TNI that was introduced in Ref. 
llOl. We summarize the strength constants and other numbers employed in the present investigation in Table [I] The 
cutoff A is the only parameter in the TNI and it is set to A = 4.7 by fitting the 3 H binding energy using the AV8' + 
TM'(99) interactions. We add the Coulomb potential to the Hamiltonian perturbatively, that is only in the P-space. 

The NCSM calculations depend in general on the HO frequency and the basis size. The basis size is limited by 
the available computational resources and the algorithm development. In this paper, we were able to reach the 6hQ 
basis space for A = 6, 7. For higher masses, we limited ourselves to the AMI basis space. The optimal HO frequency 
is selected in the NCSM by the requirement of the least dependence of the eigenergies on the frequency. Typically, 
for the p-shell nuclei we select the optimal frequency by finding the ground-state energy minimum in the largest 
accessible basis space. Due to the complexity of the calculations with the TNI, in this study we rely on our previous 
investigations done without the TNI to select the optimal frequency and perform comparative calculations with and 
without the TNI. More detailed frequency dependence investigations will be done later. 

The steps (i) and (ii) of the calculations as described in Section^ are performed using the Jacobi-coordinate many- 
body effective interaction code (manyeff) extended to handle the TNI as described in Ref. When calculating 
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the three-body effective interaction induced by the two-nucleon potentiai, V^g 123 Q, we are able reach large three- 

(3) 

nucleon basis spaces, e.g. A'max = 40, to obtain high precision three-nucleon solutions. Due to the complexity of 
the calculations with the TNI we were, however, limited in this work to smaller spaces, up to -/Vmax = 30 |5j, when 
the TNI was included in the Hamiltonian ©. This was sufficient to obtain the 3 H binding energy within 50 keV of 
the exact result, see Fig. 2 of Ref. Nevertheless, it would still be useful to further increase iVmax in the future. 
In particular because here we employ a lower HO frequency for which we expect a slower convergence than for that 
used in the 3 H calculations. In order to calculate the V^^ 23 © it is essential to use the same N^lx to obtain © 
and J7J. Therefore, we performed three calculations, one up to N^l* = 40 to calculate V^| [ 123 (|7J) in all relevant 

three-nucleon channels and two for ijfJJl and Q) in the same basis space up to = 30 for the three-nucleon channels 

up to J = | to calculate V^g^g ©. We will improve on this in our future studies as also discussed in the concluding 
section. The step (hi) is performed by a specialized code that reads the outputs of the manyeff code and the step (iv) 
is then performed using either the Many-Fermion Dynamics shell model code |16( generalized to handle three-body 
interactions or using the newly developed code REDSTICK [l7| . 

In order to get the most accurate assessment of the effect of the genuine three-nucleon interaction in the following 
we always compare calculation with and without the TNI using identical basis size, HO frequency and the clustering 
approximation, i.e., the three-body effective interaction. 

We note that bare nucleon charges and unrenormalized transition operators were used in all the transition calcula- 
tions presented in this paper. For some transitions, E2 in particular, the operator renormalization could be important. 
The operator renormalization can be calculated as outlined, e.g. in Ref. and is now under development. In the 
Gamow- Teller transition calculations, the §a constant was not included in the presented B(GT) values. 

Finally, let us remark that all the presented calculated states in this paper are the p-shell or the 0MI dominated 
states. The opposite parity and the 2?iSl-dominated states were investigated within the NCSM, e.g. in Refs. 

0, 0, |27], • Such states are not obtained at low enough energy in the basis spaces used in the present investigation. 
Therefore, in this paper comparison is made with experimental levels known to be predominantly of p-shell character. 

A. 6 Li, 6 He 

We investigated A = 6 nuclei using two-nucleon interactions in Ref. |2(i| and 6 Li in particular in Ref. Q where we 
used three-body effective interaction derived from the AV8'. Our 6 Li results with the TNI are summarized in Table ITT1 
and Fig. ^ We show the basis size dependence up to 6hil using the optimal frequency Tifl = 14 MeV. It is apparent 
from Fig. ^that the basis size dependence is quite similar in the calculations with and without the TNI and it is 
important that the excitation energies of the low-lying states in particular show a good stability and a convergence 
pattern when the 6hQ basis space is reached. At the same time we can see a significant sensitivity of the excitation 
spectra to the presence of the genuine three-nucleon interaction. It should be noted that the 3+0 state excitation 
energy decreases significantly and it gets closer to experiment when the TNI is included. Also, the 3+0 <-> 2 + 
splitting, which is underestimated with just a two-nucleon interaction, increases significantly. Our calculated 6 Li 
binding energy increases by more than 2.5 MeV when the TNI is included. With tighter binding with the TNI the 
point-proton radius decreases slightly similarly as the B(E2) value of the if — * 3 + transition. The ground-state 
quadrupole moment and the Ml transitions are not much affected, although there appears to be an improvement for 
the 2 + l -> 1^0 Ml transition. 

Our 6 He results are summarized also in Tabled Apart from increased binding energy, there are little differences 
between the calculations with and without the TNI. In the same table, we show the B(GT) values obtained for the 
6 He— > 6 Li ground-state to ground-state transition. Clearly, the TNI shifts the result in the right direction closer to 
experiment although a discrepancy still remains. Recently, Schiavilla and Wiringa found that the AV18/Urbana-IX 
interaction also over-predicts the 6 He^ 6 Li B(GT) value j24j. In fact, their obtained Gamow- Teller matrix elements 
(= y / B(GT)), 2.254(5) and 2.246(10) using two types of wave functions, compare well with our result 2.283 for 
AV8'+TM'(99). Our corresponding AV8' result is 2.305 which is basically identical to our CD-Bonn result presented 
in Ref. HJ. 

B. 7 Li, 7 Be 

Our 7 Li results with the TNI are summarized in Table ITTT1 and Fig. [21 We show the basis size dependence up to 
6hft using the optimal frequency Ml = 14 MeV. We note that the 6h£l 7 Li calculation is the most complex one among 
those presented in this paper. The three-body effective interaction basis space is defined by A% max = 9 and there are 
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741 823 056 three-nucleon M-scheme matrix elements that we input in the shell model code. Although the dimension 
of the Hamiltonian matrix is only 663 527, it is significantly less sparse when a three-body interaction is present, e.g., 
there are 4 555 058 857 non-zero matrix elements compared to 235 927 305 for a calculations with just a two-body 
interaction. 

From our 7 Li results we can draw similar conclusions as from the 6 Li results. A particularly important observation 
is the excitation energy convergence pattern, i.e., good stability of the spectra when increasing the basis size from AMI 
to QMl for both calculations with the TNI and without the TNI. There is an expected increase of binding energy and 
sensitivity of the spectra to the presence of the TNI. The level ordering is unchanged for the 9 levels that we show. 
Note that when compared to the old experimental evaluation of Ref. we obtain the reversed order of |„ | and 
| 2 i states, consistently with previous NCSM studies, e.g., (2(|. However, the new evaluation [22] changes the order 

and energies of the two levels and, in addition, adds a new | \ level. Our calculated level ordering then matches the 
new evaluated level ordering. Note that our calculation with the TNI shows a slightly improved relative level spacing 
for those three levels compared to experiment. Similarly, the relative level spacing of the | x \, |j \ and \ states 
is improved compared to experiment in the calculation with the TNI. Also we note that the lowest two excited states, 
| i and move closer to experiment when the TNI is present. 

In Table IIIII we also show our results for 7 Be and, in particular, the Gamow- Teller transitions from the 7 Be 
ground state to the ground state and the first excited state of 7 Li. There is a slight overall improvement when the 
TNI is included. The same transitions were calculated by Schiavilla and Wiringa using the the AV18/Urbana-IX 
interaction |24| . Our Gamow- Teller matrix element results compare well with theirs in particular for the ground-state 
to ground-state transition. 

C. 8 Be 

We investigated 8 Be with an emphasis on the intruder states using two-nucleon interactions in Ref. |27| and also 
in Ref. 0. Our 8 Be results with the TNI are summarized in Table llVl and Fig. [3] We compare results obtained 
using three-body effective interaction derived from the AV8' two-nucleon interaction and from the AV8'+TM 7 (99) 
interaction, respectively, in the Ahft basis space using the optimal frequency htt = 14 MeV. It should be pointed out 
that our excitation spectra calculation without the TNI is actually in a very reasonable agreement with experiment. 
It is comforting that turning on the TM'(99) three-nucleon interaction does not really change the excitation spectrum 
very much and the quality of agreement with experiment is about the same. With the TNI we gain almost 4 MeV 
of binding energy. Also, we observe a stronger isospin mixing of the 2 + and the 1 + T = 0, 1 doublets. We note that 
in experiment, the 2 + states at 16.6 and 16.9 MeV have strongly mixed T = and T = 1 components. In Table llVl 
we present the intruder 0^~0 state that was subject of our study in Ref. [27|. The position of this state, which is 
dominated by higher than OhQ components in the wave function, is strongly dependent on the basis size and still a 
much larger basis is needed for convergence. In Ref. [27J, we investigated this state using basis spaces up to 10HQ, 
although employing just the two-body effective interaction. We note that the TNI affects this state by increasing 
its energy from 26.5 MeV by about 1 MeV in the Ahtt space. However, this is not significant as we extrapolate the 
position of this state with the basis size increase to be in the 10-15 MeV range p7j . 

D. 10 B, 10 Be, 10 C 

Detailed A = 10 NCSM calculations using realistic two-nucleon interactions were presented in Ref. |2^|. In addition, 
we investigated 10 B using three-body effective interaction derived from the AV8' NN interaction in Ref. Q|. Our 10 B 
results obtained with the TNI are compared to the two-nucleon interaction calculation in Table and Fig. 0] We 
present only experimental states known to be predominantly of p-shell character as discussed earlier. In 10 B there are 
sd-shell dominated low-lying states, l + at 5.18 MeV and + l at 7.56 MeV that we obtain at much higher excitation 
energies in the basis spaces employed in this investigation j2^|. Note that in Fig. 0] the two-nucleon interaction 
excitation energy results are displayed relative to the 1^~0 state while in Table^the same results are given relative to 
the 3+0 state. We note that our calculations with realistic two-nucleon interactions predict an incorrect ground-state 
spin in this nucleus i.e., 1+0 instead of the experimental 3^~0. The same prediction is obtained by the GFMC 

method 0- Therefore, it is particularly interesting and important to investigate this nucleus using a Hamiltonian 
that includes a realistic three-nucleon interaction. Our results obtained using the TM'(99) TNI indeed show that 
the three-nucleon interaction has a positive impact on the excitation spectrum and corrects the level ordering by 
bringing the 3 + state below the l + state in agreement with experiment. Level ordering improves also for other 
excited states. This is particularly apparent when we present the excitation energies relative to the 3 + state as done 
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in Tabic Ivl In the GFMC calculations of Ref. the correct 10 B ground state was obtained using the Illinois TNI 
but the l + ground state remained when the Urbana IX TNI was employed. Combined with our observed sensitivity 
for higher excited states this shows that nuclear structure will play an important role in determining the form and 
the parametrization of the three-nuclcon interaction that, unlike the two-nucleon interaction, is still not well known. 
Concerning the binding energies, with the TM'(99) we gain almost 6 MeV for the 3 + state and about 4 MeV for the 
l + state. 

Our 10 Be and 10 C results are presented also in Table E] In particular, we paid attention to the Gamow- Teller 
transitions from the 10 B ground state to the 2+1 states of 10 Be and the transition from the 10 C ground state to 
the l + state of 10 B. Overall, we observe a clear improvement compared to experiment in calculations that include 
the TNI. In the excitation spectra, we note the increased splitting between the 2+1 and the 2 + l states when TNI is 
included, again an improvement compared to experiment. 



E. n B, n C 

There were no published NCSM calculations for 11 B up to now. Our n B results with and without the TNI are 
summarized in Table IVII and Fig. |3J We compare results obtained using three-body effective interaction derived 
from the AV8' two-nucleon interaction and from the AV8'+TM'(99) interaction, respectively, in the 4MI basis space 
using the same optimal frequency as for 10 B, i.e., Ml = 15 MeV. All shown calculated states are 0fr£l dominated. 
It is not straightforward to make a correct correspondence to the experimental T = \ levels in particular above 10 
MeV of excitation energy. We follow Refs. |3^, [3] m making assignments in Table IVII and Fig. [SJ The T = § 
level energies are taken from Refs. [U|3j|. Similarly as for 10 B, in Fig. \5\ the two-nucleon interaction excitation 
energy results are displayed relative to the \ x \ state while in Table IVTl the same results are given relative to the 
\,\ state. Given the complicated situation in neighboring 10 B, it is quite interesting to investigate 11 B. Although we 
have not done a detailed frequency dependence study as for 10 B, our present calculation suggest that the AV8' NN 
potential by itself might produce the incorrect ground state \ x \ instead of the experimentally observed \ x \ state. 
More calculations are needed to explore this issue but based on our present calculations we can definitely expect that 
a realistic two-nucleon interaction will predict the | and \ states almost degenerate. In our AMI calculation with 
the AV8' we obtain an incorrect ground state spin for 11 B. By including the TM'(99) this is corrected together with 

the level ordering of |- ^ and | 2 \ states. In addition, the level spacing improves for majority of the low-lying states 
shown in Fig. Eland Table IvTl The binding energy gain is about 6 MeV. 

As for n B, the calculation with the TNI produces a superior excitation spectrum and the correct ground-state spin 
for 11 C. An exception where there is a shift away from experiment in the calculations with the TNI appears for the 
ground-state magnetic moments of n B and n C. Recently, the Gamow- Teller transition strength for the 11 B— > n C 
was measured from the 11 B(p, n) 11 C reaction [sij . In Table ED we summarize the experimental and our calculated 
Gamow- Teller transition results. Overall, much better agreement is obtained in the calculation with the TNI. 



Detailed 12 C NCSM calculations using realistic two-nucleon interactions were reported in Ref. [l|. Here we extend 
those calculations by including the TNI. In Table IVIII we summarize some of our results. In Fig. EJ we compare 
the excitation energies obtained with and without the TNI. In addition to an increase of binding energy by 6 MeV, 
we observe a sensitivity of the low-lying spectra to the presence of the realistic three-body force and a trend toward 
level-ordering improvement in comparison to experiment. The sensitivity is the largest for states where the spin-orbit 
interaction strength is known to play a role. Note the correct ordering of the l + «-> 4+0 states and ordering and 
spacing improvement of the lowest T=l states. In fact, the lowest T = 1 state, the isospin analog state of the 12 N and 
12 B ground state comes out wrong, i.e., 0+1 instead of 1 + 1, in the calculation with just the AV8'. Once the TM'(99) 
is included we obtain the correct 1 + 1 state as the lowest T = 1 state and, moreover, the correct ordering of the next 
two lowest T — 1 states, 2 + l and + l. 

We note that in Fig. [fj]we also show low-lying 12 C states that are dominated by higher than p-shell configurations. 
We obtain such states at much higher excitation energies in the basis spaces employed in this investigation. 

The NCSM calculated binding energy of 12 C decreases with the basis size enlargement p|. We extrapolate the 
binding energy for the AV8' NN potential to be less than 80 MeV. Consequently, we expect the binding energy for the 
AV8'+TM'(99) also decrease if we were able to reach the QMl or larger basis space. Therefore, the good agreement 
of the AV8'+TM'(99) binding energy in the ATiil space with experiment shown in Table IVTT1 could be accidental. To 
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have a better insight on this issue, we need to perform more calculations, i.e., frequency and basis size dependence. 

It is well known that the isovector Ml transition from the 12 C ground state to the 15.11 MeV 1^1 state is very 
sensitive to the strength of the spin orbit interaction. Our B(M1; + — > 1 + 1) results are presented in Fig. [7| 
The calculations with the 2-nucleon interaction show saturation and under-predict the experiment by almost a factor 
of three. By including the three-nucleon interaction, the B(M1) value increases dramatically. This improvement is 
almost entirely due to the improved strength of the spin-orbit splitting when the three-nucleon interaction is included. 
The basis size dependence is stronger in this case. Clearly, if we were able to extend our calculations to the 6hfl 
space, we would get much closer to the B(M1) experimental value. We also note a substantial increase and overall 
better agreement with experiment for the B(E2; 2+1 — ► 0^0) in the calculation with the TNI (see Table IVTT|l . 

Our results for 12 N and 12 B presented in Table IVlTl further support the above discussed observations. In particular, 
the correct ground-state spin is obtained only in calculations that include the TNI. There is a significant improvement 
in level spacing even for the higher excited state, e.g., 2^1, 1^1 and 3^1. It is interesting to note the complicated 
situation concerning the magnetic moments of the 12 N and 12 B ground states. In the calculation, there is a significant 
cancellation between orbital and spin contributions and, in particular, between proton and neutron spin contributions. 
Clearly, the calculations with the TNI improve agreement with experiment, although a significant discrepancy remains. 
The Gamow- Teller ground-state to ground-state transitions behave as the Ml isovector transition from the 12 C ground 
state to the 15.11 MeV 1^1 state. The improvement compared to experiment is quite dramatic when the TNI is 
included. 



NCSM calculations for 13 C using a two-nucleon interaction were reported in Ref. |3fij |. In that work, the NCSM 
one-body transition densities were used for the DWBA description of proton and 3 He scattering on 13 C. Our present 
13 C results obtained with and without the TNI are compared in Table llXl and Fig. 00 We used the AHSl basis space 
using the same optimal frequency as for 12 C, i.e., fift — 15 MeV. The excitation spectrum of 13 C is quite complex. 
As we discussed earlier our calculated states presented in this study are all OhQ dominated. It is important to make 
a correct correspondence to the experimental levels as it is to be expected that the sd-shell dominated states would 
appear at low energies in 13 C. The character of low- lying 13 C states was studied in Ref. [37| . We use assignments 
found in Ref. in our Table 17x1 and Fig. |H1 Despite the complexity of the 13 C excitation spectra the NCSM 
calculation obtains the lowest six states in correct order in both calculations with and without the TNI. When the 
TNI is included the lowest five states excitation energies improve compared to experiment. Note in particular an 
increase of the excitation energy of the | | state. Also, the position of the lowest T — | state |. | improves. From 
Table llXl we can see that our AV8'+TM'(99) ATiQ calculation slightly overbinds 13 C. Similar comment as in Subsection 
IIIIFI on 12 C applies. We expect the binding energy decrease with the basis enlargement. 



IV. EFFECTS OF THE THREE-NUCLEON INTERACTION 



Overall, we find that three-nucleon interactions, and in particular the TM'(99) force employed in this work, do 
have an important effect on nuclear structure. Of course, the primary expectation from previous studies of three- and 
four-nucleon systems with exact methods, is that the TNI will provide more binding energy in many-body systems. 
This requirement was certainly verified with previous NCSM and GFMC studies using realistic NN interactions. This 
is also borne out in the tables, and in the discussion that follows we focus primarily on the A < 10 nuclei, where 
the results with just the NN-interaction are reasonably converged, and are likely to be within 1 MeV of the exact 
results A]. For 6 Li, the NN-interaction accounts for nearly 89% of the binding energy, while for 7 Li, 8 Be, and 10 B 
it accounts for 84-86% of the total binding. The TM'(99) force then contributes an additional binding of the order 
8-10% of the value obtained from the NN-interaction alone. What is somewhat remarkable is that the results with 
the TNI in 6 Li come quite close to the experimental binding energy (undcrbinding by less than 1 MeV), while the 
TNI significantly underbinds the A > 7 nuclei by 3.5-4.2 MeV. It is interesting that the underbinding becomes more 
dramatic with the addition just a single nucleon, i.e., from 6 Li to 7 Li. The nature of the underbinding requires further 
research. We plan to address three avenues of investigation. The first is the effects of the smaller three-particle model 
space on the convergence of the three-body effective interaction, which may be more severe for smaller oscillator 
frequencies. Although, our expectation is that the full effect of the underbinding is not due to a convergence issue. 
Efforts to upgrade the computer programs used to evaluate V^g^ NN to use larger three-particle model spaces are 
underway. Second, we are also in the process of improving our shell-model codes to extend the calculations up to 
Amax = 8 for A = 6, 7 and A max = 6 for 8 < A < 16 nuclei. This extension is important in order to be able to fully 
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assess the effect of the TNI on the structure of 6 < A < 10 nuclei. Finally, with these new capabilities, we hope to 
turn to the form of the TNI itself, which is not yet fully determined. 

The other important feature of the TNI is the enhancement of spin-orbit effects. This is particularly evident in 10 B, 
where the inversion of the 1+0 and 3 + states is a clear indication that the NN-interaction by itself is insufficient. A 
similar ground-state level inversions are observed for 11 B, and 12 N, 12 B. It is to be noted that the correct ordering is 
obtained with a realistic TNI with no free parameters (other than a cutoff fixed by reproducing the A = 3 binding 
energies). In addition, several calculated Gamow- Teller and Ml transitions are in better agreement with experiment 
when the TNI is included. These results are promising in that the inclusion of the TNI can improve our overall picture 
of the structure of light nuclei. 

V. CONCLUSIONS 

We extended the ab initio no-core shell model to include a realistic three-body interaction in calculations for p-shell 
nuclei. We presented results of first applications to 6,7 Li, 6 He, 8 Be, 10 > 1:L > 12 B, 12 N and 10 > 12 > 13 C using the Argonne V8' 
NN potential and the chiral-symmetry-based Tucson-Melbourne TM'(99) three-nucleon interaction. In general, in 
addition to increase of binding energy, we observed a substantial sensitivity of the low-lying spectra to the presence of 
the realistic three-nucleon interaction and a trend toward level-ordering and level-spacing improvement in comparison 
to experiment. The sensitivity is the largest for states where the spin-orbit interaction strength is known to play 
a role. In particular, with the TNI we obtain the correct ground-state spin for 10 ' n B and 12 N, 12 B, contrary to 
calculations with NN potentials only. In addition, we observe an overall improvement of the Gamow- Teller transition 
description in calculations with the three-nucleon interaction. The most striking case appears in A = 12 nuclei for 
the 0+0 — » l^T transitions. 

Unlike two-nucleon interaction, a detailed form and parametrization of the TNI is not well established yet, and it 
is highly likely that the sensitivity of the excitation spectra to the presence of the three-nucleon interaction could be 
exploited in determining the structure of the three-nucleon interaction itself. This is a topic that we plan to pursue 
in the future. 

Until now calculations with a realistic TNI in the p-shell has been performed only using the GFMC method 
IE 01 an d by the coupled cluster method for 16 [3!| . Our NCSM approach has some advantages compared 
to those methods such as access to wider range of nuclei and also obtaining a more complete set of excited states. 
An important issue for the NCSM is the convergence which implies the need to use as large basis space as possible. 
With some algorithm and coding developments it should be feasible to perform 8hil calculations with the TNI for the 
light p-shell nuclei and QMl calculations for the rest of the p-shell in a near future. Another important topic is the 
use of a greater variety of three-nucleon interactions and their different parametrizations including implementation of 
additional terms 0, |4(j as well as the new effective field theory interactions developed in Ref. [41] and applied in 
Ref. 42]. A closely related issue is an improvement of efficiency and accuracy of the three-nucleon system solution 
calculations needed to construct the effective interaction as described in step (i) of Section|nJ A work in this direction 
is now under way |43| . 
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FIG. 1: Calculated positive-parity excitation spectra of Li obtained in 2HQ-6HQ basis spaces using three-body effective 
interactions derived from AV8' NN potential and AV8' NN potential plus TM'(99) three- nucleon interaction, respectively, are 
compared to experiment. The HO frequency of hQ, — 14 MeV was used. The experimental values are from Ref. [21|. 
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FIG. 2: Calculated negative-parity excitation spectra of 7 Li obtained in 2ftf2-6fif2 basis spaces using three-body effective 
interactions derived from AV8' NN potential and AV8' NN potential plus TM'(99) three- nucleon interaction, respectively, are 
compared to experiment. The HO frequency of h£l — 14 MeV was used. The experimental values are from Ref. [22|. 
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FIG. 3: Calculated positive-parity excitation spectra of Be obtained in Ah£l basis space using three-body effective interaction 
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FIG. 7: Experimental and calculated B(M1) values, in /i^j, for the + — > 1 + 1 transition in 12 C. Results obtained using the 
AV8' and AV8'+TM'(99) interactions in basis spaces up to 4hQ are compared. The HO frequency of hQ, — 15 MeV was used. 
The experimental values are from Ref. |3(J. 
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3.563 


3.886 


3.933 


-E x (2+0) 


4.312 


5.010 


4.733 


EU2+1) 

x \ 1 / 


5.366 


6.482 


6.563 


£ x (l+0) 


5.65 


7.621 


6.922 


£ x (2+l) 




10.693 


10.095 


£*(1+1) 




11.525 


10.881 


B(E2;l+0 -> 3+0) 


21.8(4.8) 


7.093 


8.412 


B(M1;0 + 1 -» 1+0) 


15.42(32) 


15.374 


15.234 


B(E2;2+0 -» 1+0) 


4.41(2.27) 


3.129 


3.491 


B(M1;2+1 -> 1+0) 


0.150(27) 


0.113 


0.027 


6 He 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 




6ftO 


6hQ 


\E gs |(0+1) 


29.269 


28.189 


25.488 


r p [fm] 


1.72(4) 


1.707 


1.743 


£x(0+l) 


0.0 


0.0 


0.0 


£ x (2+l) 


1.8 


2.598 


2.638 


£ x (2+l) 




6.892 


6.231 


S x (l+1) 




7.737 


7.035 


£ x (0+l) 




10.855 


9.467 


6 He^ 6 Li 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 




6hQ 


6hQ, 


B(GT;0+1 -+ 1+0) 


4.728(15) 


5.213 


5.314 



TABLE II: Experimental and calculated energies, in MeV, ground-state point-proton rms radii, quadrupole and magnetic 
moments, as well as E2, in e 2 fm 4 , and Ml, in n%, transitions of 6 Li and 6 He as well as the B(GT) values for the 6 He 
ground state to 6 Li ground state transition. Results obtained using three-body effective interactions derived from the AV8' 
and AV8'+TM799)interactions are presented. A HO frequency of hQ, — 14 MeV was employed. The experimental values are 
from Ref. El], El El. 
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7 Li 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 


_ 


6htl 


QhQ. 


|£ g s(f~±)i 


39.270 


35.832 


32.953 


r v [fin] 


2.27(2) 


2.071 


2.104 


Q gs [e fm 2 ] 


-4.00(6) 


-2.563 


-2.679 


Mgs [a*n] 


+3.256 


+3.025 


+3.012 




0.0 


0.0 


0.0 




0.478 


0.317 


0.229 


E *(h I) 


4.65 


5.301 


5.499 




6.60 


7.388 


6.971 


S x (| 2 |) 


7.45 


8.527 


8.732 




8.75 


10.668 


10.026 


S x (| 2 |) 


9.09 


11.145 


10.793 


£x(f 2 5) 


9.57 


11.568 


11.331 


S *(fl 1) 


11.24 


12.583 


12.613 


7 Be 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 


- 




QhQ. 


i^(ri)i 


37.600 


34.176 


31.318 


r p [fm] 


2.36(2) 


2.225 


2.262 


Q gs [e fm 2 ] 




-4.305 


-4.614 


Mgs [Mn] 


-1.398(15) 


-1.153 


-1.138 


S *(fi 1) 


0.0 


0.0 


0.0 




0.429 


0.309 


0.225 


S *(fi I) 


4.57(5) 


5.232 


5.429 


S *(f 1 3) 


6.73(10) 


7.329 


6.912 




7.21(6) 


8.333 


8.514 


7 Be-> 7 Li 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 






QhQ. 




1.300 


1.401 


1.445 


B(GT;| 1 i -s- ± x |) 


1.122 


1.214 


1.231 



TABLE III: Experimental and calculated energies, in MeV, ground-state point-proton rms radii, quadrupole and magnetic 
moments, of 7 Li and 7 Be as well as the B(GT) values for the 7 Be to 7 Li transitions. Results obtained using three-body effective 
interactions derived from the AV8' and AV8'+TM^(99) interactions are presented. A HO frequency of hil = 14 MeV was 
employed. The experimental values are from Ref. |22t I23I I2B) . 
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8 Be 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 




4htl 


4hQ, 


\E gB \ 


56.50 


52.328 


48.454 


E x (0+0) 


0.0 


0.0 


0.0 


E x (2+0) 


3.04 


3.724 


3.652 


EJitO) 


11.40 


12.713 


12.402 


£ x (2+l)° 


16.63° 


16.831° 


16.630 


E*(2+0) a 


16.92° 


17.009° 


16.998 


S x (lfl) 


17.64 


18.049 


17.649 




18.15 


17.895 


17.463 


£ x (3+l) 


19.01 


19.797 


19.608 


&(3^0) 


19.24 


20.163 


19.028 


E40+1) 




20.910 


19.525 


£ x (4+0) 


19.86 


20.942 


21.144 


£ x (0+0) 




27.705 


26.478 


£ x (0+2) 


27.49 


28.931 


28.335 



°T=0 and T=l components strongly mixed 

TABLE IV: Experimental and calculated energies, in MeV, of 8 Be. Results obtained using three-body effective interactions 
derived from the AV8' and AV8'+TM'(99) interactions are presented. A HO frequency of h£l = 14 MeV was employed. The 
experimental values are from Ref. |2l|. 
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10 g 


Exp 


AV8'+TM'(99l 


AV8' 


basis space 




AMI 


AhQ, 


\E(3 + 0)\ 


64.751 


60.567 


54.833 


r„(3 + 0) ffml 


2.30(12) 


2.168 


2.196 


n,, n \p fm 2 l 
L c 1111 J 


+8.472(56) 


+5.682 


+5.937 


ll,r>-l-ri \ //ivrl 


+1.8006 


+1.847 


+1.857 


|£(1+01 

1 \ w / 1 


64.033 


60.227 


55.979 


// , _i_ „ f / / ? T 1 

An+o LMNJ 


+0 63( , 12 , 1 


+0.802 


+0.843 


x Wi / 


0.0 


0.0 


0.0 


£L(1+01 
x V 1 / 


0.718 


0.340 


-1.146 




1.740 


1.259 


1.039 


EJltO) 


2.154 


1.216 


1.664 


-Ex(2+0) 


3.587 


2.775 


1.579 


-£ x (3+0) 

X V 2 / 


4.774 


5.971 


4.363 


-Ex(2+1) 

x V 1 "V 


5.164 


5.182 


4.553 


£x(2+0) 


5.92 


3.987 


3.470 


S x (4+0) 


6.025 


5.229 


4.732 


£x(2+l) 


7.478 


7.491 


5.741 


B(E2;l+0 -» 3+0) 


4.13(6) 


1.959 


3.568 


B(E2;l+0 -» 3+0) 


1.71(26) 


1.010 


0.047 


B(E2;l+0 -» 1+0) 


0.83(40) 


3.384 


2.311 


B(E2;3+0 -» 1+0) 


20.5(26) 


3.543 


3.289 


10 Be 


Exp 


AV8' '+TM' 7991 


AV8' 


basis space 






4hti 


l-EU 


64.977 


61.387 


55.840 


r p [fm] 


2.24(8) 


2.087 


2.113 


£x(0+l) 


0.0 


0.0 


0.0 


£x(2+l) 


3.368 


3.877 


3.463 


£ x (2+l) 


5.958 


6.241 


4.706 


^x(l^l) 




8.532 


7.582 


£ x (3+l) 




9.856 


8.190 


£x(2+l) 


9.4 


10.036 


9.040 


io B ^io Be 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 


_ 


4ftfi 


4AQ 


B(GT;3+0 -» 2+1) 


0.08(3) 


0.066 


0.062 


B(GT;3+0 -> 2+1) 


0.95(13) 


1.291 


1.554 


10 c 


AV8'+TM'(99) 


Exp 


AV8' 


basis space 




47j.fi 


4AQ 


1 6 a 1 


60.321 


56.626 


51.141 


r p [fm] 


2.31(3) 


2.246 


2.279 


S x (0+1) 


0.0 


0.0 


0.0 


£ x (2+l) 


3.354 


3.913 


3.508 


£ x (2+l) 




6.133 


4.612 


^x(l^l) 




8.437 


7.453 


£x(3+l) 




9.773 


8.145 


£ x (2+l) 




10.049 


9.009 


less 10 C-, 10 B 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 






AMI 


B(GT;0+1 -> 1+0) 


3.44 


4.331 


4.748 



TABLE V: Experimental and calculated energies, in MeV, ground-state point-proton rms radii, the quadrupole moments, as 
well as E2, in e 2 fm 4 , transitions of 10 B, 10 Be and 10 C as well as selected Gamow- Teller transitions. Results obtained using 
three-body effective interactions derived from the AV8' and AV8'+TM'(99) interactions are presented. A HO frequency of 
hQ. = 15 MeV was employed. The experimental values are from Ref. pll l23t I2BI l2*sjl . 
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Exp 


AV8'+TM'(99) 


AV8' 


basis space 


_ 


4htl 






76.205 


73.338 


67.214 


r p (|-|) [fm] 


2.24(12) 


2.148 


2.175 


Q(|-i) [efm 2 ] 


+4.065(26) 


+2.920 


+2.674 


Miri) [mn] 


+2.689 


+2.176 


+2.708 




74.080 


72.816 


67.413 


[mn] 




-0.435 


-0.505 




0.0 


0.0 


0.0 




2.125 


0.522 


-0.198 


-^ x (fi 1) 


4.445 


3.565 


2.642 


-E x (| 2 i) 


5.020 


3.840 


2.492 


-Ex(5 1 5) 


6.743 


5.481 


4.946 


£x(f „ ^) 
^ 2 2 ^ 7 


8.92 


9.141 


8.130 


-E"x(§ 2 ^) 




9.819 


10.372 






10.643 


10.280 


•Kx(f g |) 


11.64 


11.989 


9.572 






12.839 


12.264 




12.916 


13.479 


12.443 


s *(§ 1 1) 




13.402 


11.139 






14.026 


12.143 


^(|-|) 


15.29 


16.588 


14.735 




16.50 


17.056 


15.968 


B ( E2 'li 1 I1 1) 


2.6(4) 


1.141 


0.522 


n c 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 


_ 


4hQ 


AhQ. 




73.440 


70.618 


64.515 


Q(|-i) [efm 2 ] 


+3.327(24) 


+2.363 


+1.627 


Mia) [Mn] 


-0.964(1) 


-0.460 


-0.923 




71.440 


70.093 


64.712 






+0.809 


+0.882 


-Mil I) 


0.0 


0.0 


0.0 




2.000 


0.525 


-0.197 


^x(ia) 


4.319 


3.584 


2.656 


S x (| 2 i) 


4.804 


3.852 


2.498 


I) 


6.478 


5.363 


4.848 




8.420 


8.943 


7.978 


n B-> n C 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 








D l u 1 '2l 2 2 1 2 J 


0.345 


0.315 


0.765 


B(GT;| 1 i -s- i x i) 


0.399 


0.591 


0.909 


r>yvj i, 2l 2 212/ 


0.961" 


0.517 


0.353 


TD/r~ ,r r. 3 — 1 » 3 — 1 \ 
B l l,1 '2l 2^ 22 2) 


0.961 a 


0.741 


0.531 


D v u L ' 2 1 2 22 2' 


0.444 6 


0.625 


0.197 



"Unresolved doublet, _E X = 4.32 + 4.80 MeV, approximately equal strength. 

b Unresolved doublet, £ x = 8.10 + 8.42 MeV, most strength in the 8.42-MeV transition. 

TABLE VI: Experimental and calculated energies, in MeV, ground-state point-proton rms radii, quadrupole and magnetic 
moments, E2, in e 2 fm 4 , transitions of 1X B and C as well as the B(GT) values for the n B to n C transitions. Results 
obtained using three-body effective interactions derived from the AV8' and AV8'+TM'(99) interactions are presented. A HO 
frequency of hQ — 15 MeV was employed. Only experimental states known to be dominantly of p-shell character are shown. 
The experimental values are from Ref. [2^. 13(1 l3ll l32| . 
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12 C 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 




Ann 




l-L'gS | 


92.162 


91.963 


85.945 


rv ffm 1 
V L J 


2 35f'2') 


2.191 


2.209 


n„, f P fm 2 l 
^2+ L c 1111 J 




4.288 


4.613 


S x (0+0) 


0.0 


0.0 


0.0 


£ x (2 + 0) 


4.439 


3.603 


3.427 


E x (l+0) 


12.710 


11.280 


13.926 


£*(4+0) 


14.083 


13.517 


12.272 


£x(l + l) 


15.110 


16.221 


16.364 


S x (2+1) 


16.106 


16.467 


17.712 


£ x (0+l) 


17.760 


17.116 


16.213 


B(E2;2+0 -> 0+0) 


7.59(42) 


4.146 


4.765 


B(M1;1+1 -> 0+0) 


0.951(20) 


0.645 


0.305 


B(E2;2+1 -> 0+0) 


0.65(13) 


0.430 


0.247 



TABLE VII: Experimental and calculated energies, in MeV, the 2^~-state quadrupole moments, as well as E2, in e 2 fm 4 , and Ml, 
in n%, transitions of 12 C. Results obtained using three-body effective interactions derived from the AV8' and AV8'+TM^99) 
interactions are presented. A HO frequency of TiQ = 15 MeV was employed. The experimental values are from Ref. [23ll30| . 
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ia B 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 


- 






1^(1+1)1 


79.578 


78.379 


72.209 


Q 1 + 1 [efm 2 ] 


1.34(14) 


+0.692 


+0.772 


Mj+i [Mn] 

L_- 


+1.003 


+0.2917 


-0.130 




0.0 


0.0 


0.0 




0.953 


0.216 


1.335 


£x(0+l) 


2.723 


0.854 


-0.158 


E K (2+1) 


3.759(6) 


3.124 


1.533 


£ x (l+l) 


5.00(20) 


4.381 


2.727 


E K (3+1) 


5.612(8) 


5.205 


3.630 


12 C ^12 B 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 


- 


4htl 


4htt 


B(GT;0 + -> 1 + 1) 


0.990(2) 


0.666 


0.255 


12 N 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 


- 


4htl 




|£(1+1)| 


74.041 


72.601 


66.466 


Q 1 + 1 [efm 2 ] 


+1.03(7) 


+0.390 


+0.386 


Mi+i 


+0.457 


+1.124 


+1.513 


S x (l+1) 


0.0 


0.0 


0.0 


£ x (2+l) 


0.960(12) 


0.236 


1.338 


£ x (0+l) 


2.439(9) 


0.825 


-0.167 


E x (2+1) 




3.104 


1.534 


E x (l+1) 




4.392 


2.733 


£ x (3+l) 




5.077 


3.547 


12 C^ 12 N 


Exp 


AV8'+TM'(99) 


AV8' 


basis space 




4HQ. 


47ifi 


B(GT;0 + -> 1+1) 


0.877(3) 


0.662 


0.254 



TABLE VIII: Experimental and calculated energies, in MeV, quadrupole and magnetic moments of 12 B and 12 N as well as 
the B(GT) values for the 12 C ground state to 12 B and 12 N ground state transition. Results obtained using three-body effective 
interactions derived from the AV8' and AV8'+TlVf799) interactions are presented. A HO frequency of HQ = 15 MeV was 
employed. The experimental values are from Ref. [251 130| . 
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13 c 




Exp 


AV8'+TM'(99) 


AV8' 


basis space 






4hQ, 


4hfl 


\E gs \ 




97.108 


99.220 


91.889 


r p [fin] 




2.29(3) 


2.170 


2.190 


Mgs \pn] 




+0.702 


+0.630 


+0.911 


-2? x (ii |) 




0.0 


0.0 


0.0 


^ x (ll 3) 




3.685 


2.501 


0.886 


-^ x (f 1 2") 




7.547 


5.781 


4.628 


2 |) 




8.860 


8.702 


8.135 






11.748 


10.054 


9.588 






12.438 


14.244 


9.795 








14.856 


12.151 


^(ll I ) 




15.108 


15.210 


13.592 


Bx (f 3 5) 






15.396 


11.441 


^ x (22 2) 






15.764 


12.756 


5) 






18.387 


14.207 


B ( E2 ;fi 2 li 




6.4(15) 


2.815 


4.224 


B(M1;|-I^i- 


1) 


0.70(7) 


0.806 


1.213 



TABLE IX: Experimental and calculated energies, in MeV, ground-state point-proton rms radii, magnetic moments, as well 
as E2, in e 2 fin 4 , and Ml, in fi%, transitions of 13 C. Results obtained using three-body effective interactions derived from the 
AV8' and AV8'+TM'(99) interactions are presented. A HO frequency of TiQ. — 15 MeV was employed. On ly e xperimental 
states known to be dominantly of p-shell character pj7| are shown. The experimental values are from Ref. [23I |3^|. 



